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\ Abstract 

We demonstrate the importance of anisotropic neutrino radiation in the mechanism of core- 
' collapse supernova explosions. Through a new parameter study with a fixed radiation field of 

neutrinos, we show that global anisotropy of the neutrino radiation is the most effective mechanism 
' of increasing the explosion energy when the total neutrino luminosity is given. We discuss the 

reason why, and demonstrate how sensitively the success of a supernova explosion depends on the 

neutrino luminosity. 

<n : 

; 1 Introduction 

> ; 

ON I For many years since the first work of Colgate & White (1966) pQ, numerical simulations of core- 
collapse supernova explosions have been exciting topics. Here a "successful" explosion must explain 
{Sj . the following observed facts: 

C*") . • The explosion energy (observed in SN 1987A): 1.5 ±0.5 x 10 51 ergs. 

O ■ 

• Asymmetry observed by polarization measurements in SN 1987A, SN 1993 J, and other super- 
CLf novae. 

2 ' • Neutron star mass: 1.4 M & . 

• 

03 . • Explosive nucleosynthesis: 56 Ni (~ 0.07 M for SN 1987A). 

>'• 

So far, no "successful" explosions have been reported in 1-D and Multi-Dimensional simulations 
except work by Wilson and Mayle (e.g. [HE]; see ) a review by Sato of this volume for the status of 
03 ■ numerical simulation studies). 

In the so-called "delayed explosion" scenario, an important role of convection, which is simulated 
with spherical neutrino radiation field, has been widely accepted. Here we study the supernova 
explosion energy with paying special attention on anisotropic neutrino emission and locally intense 
neutrino heating thereby produced. Such a new idea that anisotropic neutrino radiation may affect 
the explosion mechanism itself was first proposed by Shimizu et al. (1994) J5 a , one of the coauthors 
in the present study, and the effect of anisotropic neutrino radiation on the explosion energy was 
intensively studied by Shimizu et al. (2001) [3]. 

Figure 1 depicts the dimensionless entropy contour map |S] to explain the overall feature of this 
idea. First, we assume that a stronger neutrino flux is emitted in the direction of the rotation (polar) 
axis of a proto-neutron star than that in the equatorial plane. This anisotropic neutrino emission, for 
example, may be caused by the rotation of the core of a progenitor star (see jlj for other possibilities). 
Then high-entropy hot bubble is created along the polar axis behind the primary shock. Convective 
motion occurs, but the convection becomes very global compared with other models with spherical 
neutrino radiation. The final important feature of our model is that the shock wave itself is deformed 
as seen in Figure 1. It was emphasized in Shimizu et al. 14 that a reasonable degree of anisotropy, 
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Figure 1: Dimensionless entropy contour map jS] with schematic notations. Entropy is calculated 
for the model of l z /l x = 1.1 and T u = 4.7 MeV at t = 160 ms after the shock stall; see section 2 for 
explanation of this model. 

lz/lx = ~105 - ~1.2, is enough to increase the explosion energy by a large factor. In the above, \ z 
and l x are, respectively, the local neutrino flux in the polar axis and that in the equatorial plane, (see 
[UEl for details of parameterization of neutrino flux distributions). Shimizu et al. (2001) concluded 
that global anisotropic neutrino radiation is an alternative mechanism to revive the shock wave and 
lead to a successful explosion other than so far suggested "convective trigger" . 

In this article we are going to show that global anisotropy of the neutrino radiation is the most 
effective mechanism of increasing the explosion energy than the models of spherically symmetric and 
fluctuated neutrino radiation when the total neutrino luminosity is given. We discuss the reason 
why, and demonstrate how sensitively the success of a supernova explosion depends on the neutrino 
luminosity. 

Our simulation is performed by solving 2-dimensional hydrodynamic equations in spherical co- 
ordinates. A generalized Roe's method is employed to solve the hydro dynamical equations with 
general equations of state (EOSs). Since our method is very rigorous to treat neutrinos above the 
neutrinosphere, we put our initial condition (t = 0) when a spherical steady shock is stalled above 
the neutrinosphere. The details of our numerical technique, together with the EOS and the initial 
condition used, are described in the previous article [3]. It is noted here that we have improved the 
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Figure 2: Evolution of thermal and kinetic energy (E t h+Ek in ), gravitational energy (E grav ), and 
explosion energy (E exp i) for the model of \ z /\ x = 1.1 and T u = 4.7 MeV. 

numerical code in the papers [I] and 6 to avoid a numerical error near the pole; although the error 
was not serious for the investigation of the explosion energy, this may affect the results of explosive 
nucleosynthesis. 

2 Reasoning of Powerful Explosion with Anisotropic Neutrino Ra- 
diation 

Figure 2 shows the evolution of the explosion energy, as well as the thermal, kinetic, and gravitational 
energies for models of l z /l x = 1-1 and T u = 4.7 MeV. The quantity T„ is the neutrino temperature, 
and we consider here only electro-type neutrinos and the total neutrino luminosity is fixed (see |Hj for 
details). In Figure 2, we observe that the difference in the evolution of the explosion energy between 
the globally anisotropic model and the spherically symmetric model is clearly prominent. 

Two dashed lines in the explosion energy (E exp i) in Figure 2 are for the cases when the neutrino 
field is fluctuated in space (see 6 for details). Such neutrino fluctuations may occur due to gravita- 
tional oscillation of proto-neutron stars as proposed by Burrows et al. [Zj. The model for the lower 
dashed line of E exp i in Figure 2 has larger mode number of fluctuations in the neutrino flux. It is 
therefore seen that the explosion energy decreases as the mode number of fluctuations in the neutrino 
flux increases, and finally approaches that of spherical explosion. 

Let us mention here the reasoning of this powerful explosion for the globally anisotropic neutrino 
radiation model. Figure 3 shows the accumulation of the absorbed and emitted energies due to 
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Figure 3: Accumulation of the absorbed and emitted energies due to neutrino heating (E a b s ) and 
cooling (E em it), respectively, for the anisotropic and spherical models of T u = 4.7 MeV. Note the 
difference in the emitted energy. 

neutrino heating and cooling for both globally anisotropic = 1.1) and spherical models. 

The balance between heating and cooling determines the evolution of the explosion energy. The 
crucial point here is the following: The neutrino heating is caused by neutrino absorption due to 
nuclei and neutrino scattering off electrons and positrons, and hence the heating rate depends on 
the neutrino temperature T u : It varies as T^. On the other hand, the neutrino cooling is caused 
by neutrino emission due to electron captures and the emission of thermal (photo, pair, plasma) 
neutrinos. The cooling rate is then a very sensitive function of the matter temperature T, roughly 
proportional to T 6 . 

The reasoning of the powerful explosion for the globally anisotropic neutrino radiation model can 
be explained as follows. The locally intense neutrino radiation in our model causes local heating 
along the polar axis. The thermal pressure around the heated matter becomes high, and hence it 
pushes part of the shock wave outward. Namely, the shock revival occurs in one direction first along 
the axis of rotation. Then the pressure gradient does work along the shock front, and the shock 
revival prevails to all directions (Note the Rankine-Hugoniot condition). This leads to the earlier 
shock revival than the spherical model. Thus the matter temperature behind the shock decreases 
rapidly due to expansion. Accordingly, the neutrino cooling rate rapidly drops. Heating rate, on 
the other hand, remains almost unchanged between the the anisotropic and the spherical models. 
Heating dominates cooling as a result (see Figure 3). Finally, the explosion energy increases because 
of the suppression of neutrino cooling. 

Notice here that the position where the cooling rate suddenly starts to decrease in Figure 3 (i.e., 
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Figure 4: Same as Figure 2 except T„ = 4.65 MeV. 

t ~ 80 ms for the anisotropic, and t ~ 230 ms for the spherical model) corresponds to the shock 
revival time, which can be seen in Figure 2 as the position where the explosion energy suddenly starts 
to rise up. 

3 Sensitivity of Explosions to Neutrino Luminosities 

In Figure 4 we show a model case with a reduced temperature which fails to explode when the neutrino 
field is spherically symmetric and fluctuated, but does produce an adequate explosion when the field 
is reasonably anisotropic (\ x /l x = 1-1); the difference between the model with global anisotropy and 
the other models is extremely remarkable. 

We now see that two spherically symmetric models which yield very different explosion outcomes 
in Figures 2 and 4, but these differ by only 1 % in neutrino temperature. This is a clear demonstration 
of the sensitivity of the supernova problem on the neutrino luminosity and energy. 

4 Concluding Remarks 

It has been discussed that the globally anisotropic (jet-like) explosion is the most effective mecha- 
nism to increase the explosion energy than spherical and neutrino-fluctuated explosions when the 
total neutrino luminosity is given. While the enhancement of total neutrino luminosity is the most 
important ingredient for driving to an energetic explosion as concluded by Janka & Miiller [HIE]) 
the total neutrino luminosity cannot be simply increased to explain the observed explosion energy. 
This is because such treatment often leads to the problem of Ni overproduction, especially in the case 



of essentially spherical models. We have demonstrated that the enhancement of the total neutrino 
luminosity is not always necessary: the less luminous the neutrino radiation is, the more important 
the neutrino anisotropy is (compare Figures 2 and 4) . This is a confirmation of the claim by Shimizu 
et al. [I], using a different neutrino flux distributions and an updated numerical code. 

The degree of anisotropy assumed here, l z /l x = 1-1, can be translated into the rotation period of 
a proto- neutron star of ~ 15 ms on the assumption of the Maclaurin spheroid (see |3|). We note that 
the following saturation properties of the neutrino anisotropy has been discussed in the papers |HI§|; 
the effect of anisotropic neutrino radiation on the explosion energy appears with a small value (e.g., 
lz/L ~ 1.05), and the effect saturates with a certain value (e.g., l z /l x ~ 1-2; the actual value depends 
on the assumed form of the local neutrino flux) . We emphasize here again that a large anisotropy in 
the neutrino radiation is NOT required for the explosion mechanism. 

Finally, we also showed that the success of an explosion depends strongly on a subtle difference 
in the neutrino energy and luminosities (compare the results of Figures 2 and 4; the temperature 
difference between 4.70 and 4.65 MeV is only 1 %). What is most important is that anisotropic 
neutrino radiation, or local neutrino heating, has been suggested as an alternative key to the explosion 
mechanism of super novae. 
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